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EFFECTS OF AIRCRAFT RELATIVE DENSITY ON SPIN AND 


RECOVERY CHARACTERISTICS OF SOME 

CURRENT CONFIGURATIONS 

By William D. Grantham and Sue B. Grafton 
Langley Research Center 


SUMMARY 


An analytical study has been conducted on a high-speed digital computer 
utilizing six-degree-of-f reedom equations of motion to examine the effects of 
relative density on the spin and recovery characteristics of four configura- 
tions which are representative of modern airplanes. Two approaches were used: 
Computations were made simulating conditions for which the airplane obtained a 
disturbance that put it at a high angle of attack with applied rotation in a 
near- developed spin condition at various initial altitudes to determine whether 
a developed spin would ensue. After it was determined that developed spins did 
ensue in the first group of calculations, separate calculations were made to 
determine whether a spin could be entered starting at or near trimmed gliding 
flight (l g stall maneuvers). 

The results indicate that the effects of changing relative density on 
developed spins and recoveries were as follows: An increase in relative density 

gave faster rotating spins, higher rates of descent, lower values of the spin 
coefficient, little change in angles of attack and sideslip, and recoveries, if 
obtained, were slower. Changes in relative density can make the difference 
between a spin and a no- spin when entry is attempted by means of a 1 g stall 
maneuver, but the effect of relative density is not consistent. About the only 
generalization that can be made on the effect of relative density on the spin 
entry is that increases in relative density cause increased roll oscillations 
during the spin-entry motions. 


INTRODUCTION 


Reference 1 describes the results of model tests made more than 20 years 
ago in the Langley 20-foot free- spinning tunnel to determine the effects of 
relative density on spin and recovery characteristics of airplanes. In that 
study, the range of relative-density parameters used was 6 to 12. Over the 
ensuing years, airplane configurations have changed considerably and the range 
of values of relative density has greatly increased. An analytical study has, 
therefore, been made to examine the effects of relative density on spinning for 



four more recent configurations, and the results are reported herein. In this 
study the range of relative-density parameters covered was from 17 to ^ 87 . The 
present paper also includes consideration of spin-entry characteristics, whereas 
reference 1 was concerned only with fully developed spins and recoveries 
therefrom . 

The four configurations used in this study were considered to be representa- 
tive of modern aircraft and were as follows: a stub-wing research vehicle, a 

sweptback wing fighter, a delta-wing bomber, and a delta-wing fighter. The 
relative-density parameter was considered to be varied by changing altitude. 
Included in the study were brief calculations to determine the effect of varying 
the magnitude of some of the aerodynamic parameters for various simulated 
altitudes . 


SYMBOLS 


The body system of axes is used. This system of axes, related angles, and 
positive directions of corresponding forces and moments are illustrated in 
figure 1. 
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longitudinal-force coefficient. 
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mean aerodynamic chord, ft 

longitudinal force acting along X body axis, lb 
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side force acting along Y body axis, lb 
vertical force acting along Z body axis, lb 
acceleration due to gravity, ft/sec^ 

acceleration due to gravity at sea level, J 2.17 ft/sec^ 
altitude at beginning of time increment, ft 
altitude at end of time increment, ft 

moments of inertia about X, Y, and Z body axes, respectively, 
slug-ft 2 

rolling moment acting about X body axis, ft-lb 
pitching moment acting about Y body axis, ft-lb 
yawing moment acting about Z body axis, ft-lb 
mass of airplane, W/g, slugs 

components of resultant angular velocity about X, Y, and Z body 
axes, respectively, radians/sec 

radius of earth, 3 , 956.67 miles 

wing-surface area, sq ft 

time, sec 

components of resultant velocity V R along X, Y, and Z body axes, 
respectively, ft/sec 

vertical component of velocity of airplane center of gravity (rate of 
descent), ft/sec 

resultant linear velocity, ft/sec 
weight, lb 

longitudinal, lateral, and vertical body axes of airplane, - 
respectively 

angle of attack, angle between relative wind Vr projected into 

XZ-plane of symmetry and X body axis, positive when relative wind 
comes from below XY body plane, deg 



angle of sideslip, angle between relative wind V R and projection 

of relative wind on XZ-plane, positive when relative wind comes 
from right of plane of symmetry, deg 

total deflection of left and right ailerons with respect to each 
other, positive with trailing edge of right aileron down (left 
stick), deg 

elevator deflection with respect to fuselage reference line, positive 
with trailing edge down, deg 

rudder deflection with respect to fin, positive with trailing edge 
to left, deg 

total angular movement of X body axis from horizontal plane meas- 
ured in vertical plane, positive when airplane nose is above hori- 
zontal plane, deg 

airplane relative-density parameter, — - — 

gpSb 


air density, slugs/cu ft 

resultant angular velocity, radians/sec 

angle between Y body axis and horizontal measured in vertical plane, 
positive for erect spins when right wing is downward and for 
inverted spins when left wing is downward, deg 

total angular movement of Y body axis from horizontal plane measured 
in YZ body plane, positive when clockwise as viewed from rear of 
airplane (if X body axis is vertical, 0 e is measured from a 

reference position in horizontal plane), deg 

horizontal component of total angular deflection of X body axis from 
reference position in horizontal plane, positive when clockwise as 
viewed from vertically above airplane, deg 

incremental rolling -moment coefficient due to aileron deflection, per 
deg 

incremental rolling-moment coefficient due to rudder deflection, per 
deg 

incremental yawing-moment coefficient due to aileron deflection, per 
deg 

incremental yawing-moment coefficient due to rudder deflection, per 
deg 



incremental side-force coefficient due to aileron deflection, per deg 




incremental side-force coefficient due to rudder deflection, per deg 
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A dot over a symbol represents a derivative with respect to time. 


PROCEDURES AND CALCULATIONS 


Spin entry and developed spin motions were calculated by a high-speed 
digital computer which solved the equations of motion and associated formulas 
listed in the appendix. The equations of motion are Euler’s equations repre- 
senting six degrees of freedom along and about the airplane body system of axes. 
(See fig. 1 for illustration of body axes.) The mass and dimensional character- 
istics used in the calculations are listed in table I and the planviews of the 
four configurations designated A, B, C, and D are shown in figure 2. 

In general, the aerodynamic data used were nonlinear as shown in figures 3 
to 8. The data for these plots were obtained from references 2 to 7* No meas- 
ured values of the lateral force increments resulting from deflecting the ailer- 
ons and rudder were available for configuration B. (See figs. 5 and- 6.) Inasmuch 
as previous studies have indicated that these incremental forces have little or no 
effect on spin entry and spin characteristics no attempt was made to estimate them. 
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The oscillation-type rotary derivatives presented in figures 7 and 8 were 
obtained as combination derivatives which include the effects of 0 - that is, 
is actually + sin Cn^ is actually ^C n p + C n . sin aj, Cj r 

is actually (Ci T - C^- cos a^, C nr is actually (^n r - C n . cos aj, and so 

forth- However, inasmuch as the full derivatives could not be separated into 
their component parts, it was arbitrarily decided for this study to treat the 
derivatives as though they were due solely to angular velocities about body 
axes. The rolling moment due to yawing C^ r and the yawing moment due to 

rolling parameters were set equal to zero for configurations A and B. For 

all configurations, no effects of rolling and yawing on side force were included. 


In addition, constant values of 
lows: Configuration A, 


Qm n were used for each configuration as fol- 


C m ^ = -10; configurations B and C, 


Cmq = -2; and con- 


figuration D, 


Cm = "1 • 
1 


Two approaches were used: Computations were made simulating conditions 

for which the airplane obtained a disturbance that put it at a high angle of 
attack with applied rotation in a near-developed spin condition at each alti- 
tude (that is hQ = 15,000 feet, 30,000 feet, 45,000 feet, and 60, 000 feet) to 
determine whether a developed spin would ensue. This technique simulated that 
used in the Langley 20-foot free-spinning tunnel in which small dynamic models 
are launched in near-developed spin conditions and then freely proceed to either 
developed spins or to "no-spin" dive-out or roll-over motions, depending upon 
the design and mass characteristics. After it was determined that developed 
spins did ensue in the first group of calculations, separate calculations were 
made to determine whether a spin could be entered starting at or near trimmed 
gliding flight (l g stall maneuver). 

Spin recovery attempts were made by deflecting the rudder against the 
direction of yaw and the ailerons with the direction of yaw (left rudder and 
right stick when spinning to the pilot's right), because these control deflec- 
tions are the optimum for recovery from developed spins for airplanes loaded 
heavily along the fuselage (see ref. 8), as are the subject configurations. The 
elevators were left in the initial up position for all cases. During the present 
study, a spin is considered to be terminated when either the spin-rotation ceases 
or the angle of attack becomes and remains less than the stall angle within 10 
turns after the recovery controls are applied. Usually when the angle of attack 
becomes less than the stall angle, the airplane enters a steep dive without 
significant rotation (r « 0). In some cases, however, the airplane may be 
turning or rolling in a spiral glide or an aileron roll. Also, sometimes the 
airplane may roll or pitch to an inverted attitude from the erect spin and may 
still have some rotation, but it is out of the original erect spin. 
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RESULTS AND DISCUSSION 


The calculated results are presented in figures 9 to 1 6 as time histories 
of angle of attack a, angle of pitch 0 e , angle of sideslip p, angle of roll 
yawing velocity r, control- surface positions, and spin turns completed. 

It should be noted that the scales on these figures are not consistent and the 
reader should be careful when trying to compare various time histories. Also, 
it may be noted that attempted recoveries are indicated by dashed lines, spins 
by solid lines. 


Simulated Launching With Rotary Motion 

Elevator up, rudder right, and stick left conditions were used in each cal- 
culation to promote a spin to the pilot’s right. The initial flight conditions 
used at each relative density investigated are presented in table II. 

The calculated results for all configurations showed that developed spins 
ensued from launchings at all four altitudes; the magnitudes of the pertinent 
parameters are presented in table III. 

Configuration A .- Two representative time histories for configuration A 
are presented as figure 9* Figure 9(a) represents the motion obtained starting 
at 15,000 feet (jj, = 75) and indicates that after approximately 10 turns a 
developed spin condition had been achieved. The angle of attack was oscil- 
lating from approximately 84° to 89°, the angle of sideslip was oscillating 
from 1° to -8°, and the rate of yaw was approximately 2.8 radians per second. 
This is the point at which the recovery attempt was made, and as seen in fig- 
ure 9( a ) no recovery was achieved within 10 additional turns. Figure 9("b) 
represents the spin obtained by starting at 60,000 feet (\i = 487) and shows 
that after approximately 10 turns a ~ 85° to 88°, |3 « 0° to -4°, and 

r « 3*1 radians per second. A recovery was attempted at that point and, as 
shown, none was achieved within 10 additional turns. 

The results obtained for this configuration indicate that an increase in 
the relative density gave slightly faster rotating spins, higher rates of 
descent, and lower values of the spin coefficient £2b/2V R . The small changes 

in angle of attack and sideslip are negligible. 

Configuration B .- Two representative time histories of the spin and recov- 
ery motions obtained for configuration B are presented as figure 10. 

Pertinent results from these calculations indicate that when this configu- 
ration was launched into a near-spin condition, an increase in the relative 
density gave faster rotating spins, higher rates of descent, lower values of 
the spin coefficient, and slower recoveries. The effect of relative density 
variations on the angles of attack and sideslip experienced during these four 
spins was again considered to be negligible. (See table III.) 
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Configuration C .- Spins were obtained at all four altitudes for configu- 
ration C, from which no recoveries were achieved. Two representative time 
histories are presented as figure 11. 

The results (table III) indicate that, as in the case of configurations A 
and B already discussed, when this configuration was launched into a near-spin 
condition, an increase in the relative density gave slightly faster rotating 
spins, higher rates of descent, and lower values of the spin coefficient; how- 
ever, there was little change in the angles of attack and sideslip. 

Configuration D .- The calculated time histories for configuration D show 
that recoveries were achieved from each spin. Two representative time histories 
of the spin and recovery are presented as figure 12. 

The results shown in table III indicate that, as in the case of the con- 
figurations already discussed, when this configuration was launched into a near- 
spin condition, sin increase in the relative density gave faster rotating spins, 
higher rates of descent, lower values of the spin coefficient, and slower 
recoveries. The small changes in angles of attack and sideslip were again con- 
sidered to be negligible. 


Summation of results for launchings with rotary motion.- It was concluded 
from the results discussed thus far that developed spins ensued for all con- 
figurations and initial conditions studied, and that the effects of variations 
in relative density were similar on all of the configurations; that is, an 
increase in relative density gave faster rotating spins, higher rates of descent, 
lower values of the spin coefficient, little change in angle of attack and side- 
slip, and recoveries, if obtained, were slower. These results are considered 
to be in general agreement with those of reference 1; except that in the former 
study, somewhat larger effects of relative density on angle of attack and side- 
slip appeared to be indicated. 


Simulated Spin Entry Motion 

After it was found that, once achieved, developed spins could be maintained 
for all Gases investigated, attempts were made for each configuration to enter 
the spin starting at or near trimmed gliding flight and flying the airplane up 
through the stall angle of attack (l g stall maneuver). Back stick was used to 
stall the craft, right rudder was applied at or just after the stall to yaw the 
craft to the right, and left stick was applied variously timed, to attempt to 
promote spin entries to the right in all calculations. Calculations were made 
that simulated initial altitudes of 15 , 000 , 50 , 000 , 45,000 and 60, 000 feet for 
each of the four configurations, as was done for the computations where each 
craft was launched with applied rotation. Initial conditions used are shown in 
the second part of table II. 

Configuration A .- The calculated time histories for configuration A are 
presented as figure 15 , and, as can be seen, spins were obtained at all four 
altitudes . 
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Comparison of the spins at four altitudes (fig- 13) indicates that the only 
noticeable effect of change in relative density on the spin entry characteris- 
tics is that as the relative density is increased the roll angles experienced 
during the initial phases of the spin entry are larger and that the ensuing 
rolling and pitching motions were more oscillatory- After approximately five 
spinning turns, the spin motions achieved at the various altitudes are similar. 

Recoveries were attempted from each of these spins after seven spinning 
turns had been completed from initial Ig stall entry and, as can be seen from 
the respective time histories, it took approximately 6 to 7 additional turns to 
achieve recovery in each instance. It is believed that the reason recoveries, 
even though poor, were possible from these spins whereas, as previously dis- 
cussed, they were not obtained from the spins which ensued after this configura- 
tion was launched in a near- spin condition is that in the second group of cal- 
culations, the recovery controls were applied before the spin- rotation rates had 
increased to the magnitudes experienced in the first group of calculations. The 
lower rotation rates enabled recoveries, even though poor, to be achieved. 

Configuration B .: The computed time histories for configuration B are pre- 

sented in figure 14. The characteristics of motions obtained in the attempted 
spin entries were very different for the different altitudes, but there was no 
consistent trend in the results. For example, when a spin entry was attempted 
at an altitude of 15,000 feet (fig. 14(a)), the craft rolled over to the right 
(0e > 360 o ) and then continued to oscillate in roll to such an extent that no 
spin was obtained; whereas, for an initial altitude of 30,000 feet (fig. 14(b)), 
the craft rolled over to the right, maintained an angle of attack above the 
stall angle, began to oscillate approximately ± 15 ° in roll, and continued to 
spin. The time history shown in figure 14(c), which was computed by using an 
initial altitude of 45,000 feet, shows that the craft rolled over to the right 
twice (0 e > 700 °) and then continued to roll until such time as the angle of 
attack became less than zero. The direction of turning was steadily changing 
and therefore could not be called a spin. When a spin entry was attempted by 
using an initial altitude of 60,000 feet (fig. 14(d)), the craft rolled over to 
the right twice and began to oscillate from approximately - 50 ° to 70 ° in roll 
angle and made three spinning turns to the right before it oscillated out of 
the spin. The angle of attack became less than zero and the rate of yaw was 
approaching zero when the computation was stopped. No recovery was attempted 
from the only spin obtained (i.e., at ho = 30,000 feet). 

Configuration C .- The time histories for configuration C are presented in 
figure 15 and, as can be seen, spins were entered at each altitude. Again, as 
in the spins entered with configuration A, the computed rolling and pitching 
motions were more oscillatory for the higher values of relative density. 

Recoveries were attempted from these spins after seven turns had been com- 
pleted with the following results: (l) no recovery was achieved from the spin 

at 15,000 feet, (2) recovery was achieved from the spin at 30,000 feet in 
approximately two additional turns, ( 3 ) less than one additional turn was 
required to achieve recovery from the spin at 45,000 feet. (See figs. 15 (a) 
to 15 (c), respectively.) No recovery was attempted from the spin obtained at 
an initial altitude of 60, 000 feet. Apparently, the primary reason that 
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recoveries were obtained at the higher altitudes and not at 15,000 feet is that 
at the higher altitudes the spin motions were more oscillatory and the rotation 
rates slower. This is also the probable reason why no recoveries were achieved 
from the spin obtained by launching the craft with applied rotation (table III); 
those spins were more steady in nature and yawing at a somewhat faster rate when 
recovery controls were applied. (Compare fig. 11 with fig. I 5 .) 

Configuration D .- The calculated time histories for configuration D are 
presented in figure 1 6 . It is shown that spins were entered at each altitude 
investigated except at 60, 000 feet, wherein the craft rolled to such an extent 
that a spin entry was prevented. No recoveries were attempted for any condi- 
tion. About the only generalization that can be made about the results for 
configuration D is that as the relative density was increased the amount of roll 
experienced increased. 

Summation of results of spin ent ry ca lcul ations . - It is concluded from 
these calculations that there is no consistent effect of relative density on 
whether an airplane has more or less tendency to enter a spin. About the only 
generalization that can be drawn from the results of the spin entry calculations 
is that increasing relative density causes larger roll-angle oscillations in the 
initial phases of the spin-entry attempt. Apparently, changes in relative 
density can affect the roll characteristics to such an extent that spins may or 
may not be obtained from a 1 g stall maneuver for configurations and conditions, 
all of which would spin when launched with spinning rotation. For example, for 
two of the subject configurations, spins could be entered over the complete 
range of relative density investigated; for the third conf iguration, spins could 
be entered at all relative densities except the maximum studied. For the fourth 
configuration, spins could be entered at a medium value of relative density 
(craft rolled over to the right once, began to oscillate ± 15 ° and continued 
spinning), but when the relative density was appreciably decreased or increased, 
no spin could be entered (rolling did not stop after initial roll over). 

It appears that in order to predict the effect of changing relative density 
on spin-entry characteristics for any particular configuration, an investiga- 
tion must be made on the specific design. 


AERODYNAMIC VARIATIONS 


Inasmuch as it was found that there were some altitudes at which some con- 
figurations would not enter a spin by means of a 1 g stall maneuver, and that 
the roll-motion characteristics were very important in determining whether a 
spin could be entered, a few calculations were made to determine whether large 
increases in the values of and (arbitrary values used were twice the 

basic values) would enable spins to be entered at these altitudes. 


As mentioned previously, spins on configuration D were entered at all alti- 
tudes except the maximum (60, 000 feet). Additional calculations were made 
wherein the basic values of C 1 and C^p were arbitrarily increased by 
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various amounts in an attempt to enter a spin at a simulated altitude of 
60,000 feet with this configuration. The resulting time histories showed that 
increasing the negative values of C^ p did not enable a spin to be entered, 

whereas the larger negative values of C^p (increased effective dihedral) 

allowed a spin entry. These effects are in agreement with the results obtained 
during previous analytical studies of the spin characteristics of other delta 
wing configurations and indicated that the magnitude of Cj p had little effect 

on the initial phases of the spin entry, although the magnitude of C^p can 

affect the spin after it has developed. Also, as to the effects of C^p on 

spin entries, it has been pointed out in references 3 and 9 that the magnitude 
of C^p can determine whether an aircraft can enter a spin. 

As discussed previously, spins could be entered on configuration B at an 
altitude of 30,000 feet, but if the altitude was reduced to 15,000 feet or 
increased to ^5,000 and 60, 000 feet, no spins could be entered. Some calcula- 
tions were made wherein the basic values of C^ and C^p were arbitrarily 

increased by various amounts in an attempt to enter spins at 15, 000 feet and 
60,000 feet. The time histories computed for an initial altitude of 15,000 feet 
indicated that increasing the basic negative values of C^ p or C^p still 

would not enable a spin to be entered. However, when the increased negative 
values of both C^p and C^p were used, the craft did enter a spin. The 

results at 60, 000 feet showed that when larger negative values of C^ and/or 

C^p were used, spins could not be entered. 

These results indicate that even when relatively large negative values of 
C^p and/or C^p are used, spins cannot always be entered at any given altitude 

on all configurations. 


CONCLUSIONS 


The following conclusions are drawn from the present analytical study on 
the effects of airplane relative density on spin and recovery characteristics 
of four configurations representative of modern airplanes. 

1. Trends obtained as to the effect of changing relative density on 

developed spins and recoveries were generally similar to those noted in earlier 
experimental free- spinning tunnel model tests made for airplanes with unswept 
wings and a much smaller relative density range, and were as follows: An 

increase in relative density gave faster rotating spins, higher rates of descent, 
lower values of the spin coefficient, little change in angle of attack and side- 
slip, and recoveries, if obtained, were slower. 

2. Changes In relative density can make the difference between a spin and 
a no- spin when entry is attempted by means of a 1 g stall maneuver, but the 


11 



effect of relative density is not consistent. About the only generalization 
that can be made on the effect of relative density on the spin entry is that 
increases in relative density cause increased roll oscillations during the spin- 
entry motions . 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 3, 1964. 


APPENDIX 


EQUATIONS OF MOTION AND ASSOCIATED FORMULAS 


The equations of motion used in the calculations were: 
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In addition, the following formulas were used: 


a = tan“^ ^ 


0 = sin ^ — 
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V = -u sin 0 e + v cos 0 e sin 0 e + w cos 0 e cos 0 g 
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TABLE I.- MASS AND DIMENSIONAL CHARACTERISTICS 



16 



TABLE II.- INITIAL CONDITIONS USED IN CALCULATIONS 


Altitude, 

ft 


15. 000 

30. 000 

45.000 

60. 000 


15.000 

30. 000 

45. 000 

60.000 


15. 000 

30. 000 

45. 000 

60. 000 


15. 000 

30. 000 

45. 000 

60. 000 


Altitude, 

ft 


15,000 

30. 000 

45.000 

60. 000 


15,000 

30. 000 

45.000 

60. 000 


15. 000 

30. 000 

45. 000 

60. 000 


15. 000 

30.000 

45. 000 

60. 000 


(a) Simulated launch with rotation; p = 0 e = \jr e = v= q- O 


a. 

9 e> 

», 

v. 


P, 


deg 
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(b) Simulated spin entry; 3~0 e = '»|r e = v = p= r = O 
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TABLE III.- SOME PERTINENT RESULTS OBTAINED WEEN CONFIGURATIONS 


WERE LAUNCHED INTO A NEAR-SPIN CONDITION 
(~A11 values taken after approximately 10 turns except where noted] 
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Figure 9*- 


(a) h 0 = 15,000 ft. (b) h Q = 60,000 ft. 

Calculation simulating launching with rotary motion. Configuration A; S e = -30°. 
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(a) h 0 = 15,000 ft. 


(b) h 0 = 60,000 ft. 


Figure 10.- Calculation simulating launching with rotary motion. Configuration B; 6 e = 
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(a) h 0 = 15,000 ft. (b) h 0 = 60,000 ft. 

Figure 11.- Calculation simulating launching with rotary motion. Configuration C; 5 e = -20°. 
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Figure 12.- Calculation simulating launching with rotary motion. Configuration D 
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(a) h 0 = 15,000 ft. (b) h 0 = 30,000 ft. 

Figure 13 .- Calculated spin entry attempt. Configuration A; b e - -30°. 
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(c) h 0 = 14-5,000 ft. 


Figure 13 •- Concluded. 


(d) ho = 60,000 ft 
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Figure 14.- Calculated spin entry attempt. Configuration B; 5 e = -50°. 
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Figure 15 .- Concluded 
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figure 1 6.- Concluded. 
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